vr=evalr(epls);
for i=1:pop_s
    plo=br(pop_opt(i));
    vrp=evalr(plo);
    u(i)=abs((vr-vrp)*pop_opt(i));
end
[mu,mi]=min(u);
